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Abstract 

In this paper, we calculate the contribution of the littlest Higgs(LH) model to the 
process e^e~ — > ZHH at the future high energy e'^e~ collider (/LC). The results show 
that, within the parameter spaces preferred by the electroweak precision, the deviation of 
the total cross sections from its SM value varies from a few percent to tens percent. The 
correction of the LH model to the process might be detected at the future ILC experiments 
in the favorable parameter space. On the other hand, we find that the correction of the 
LH model is sensitive to the trilinear Higgs coupling in some case and the process can 
also provide us a chance to probe such coupling in the LH model. 
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1 Introduction 



The standard model(SM)provides an excellent effective field theory description of almost all 
particle physics experiments. But in the SM the Higgs boson mass suffers from an instability 
under radiative corrections. The naturalness argument suggests that the cutoff scale of the SM 
is not much above the electroweak scale: New physics will appear around TeV energies. Among 
the extended models beyond the SM, the little Higgs model offers a very promising solution 
to the hierarchy problem in which the Higgs boson is naturally light as a result of nonlinearly 
realized symmetry [11 [2l [3]. The key feature of this kind of models is that the Higgs boson is 
a pseudo-Goldstone boson of an approximate global symmetry which is spontaneously broken 
by a vacuum expectation value (VEV) at a scale of a few TeV and thus is naturally light. 

The most economical little Higgs model is the littlest Higgs(LH) model, which is based on 
the SU{5)/ S0{5) nonlinear sigma model jl]. It consists of a SU{5) global symmetry, which is 
spontaneously broken down to 50(5) by a vacuum condensate /. In the LH model, a set of new 
heavy gauge bosons(-BH, Zh, Wh) and a new heavy-vector-like quark(T) are introduced which 
just cancel the quadratic divergence induced by SM gauge boson loops and the top quark loop, 
respectively. The distinguishing features of this model are the existence of these new particles 
and their couplings to the light Higgs. Measurement of these couplings would verify the struc- 
ture of the cancelation of the Higgs mass quadratic divergence and prove the existence of the 
little Higgs mechanismjS]. 

The hunt for the Higgs boson and investigation of its properties is one of the most impor- 
tant goals of present and future high energy collider experiments. The precision electroweak 
measurement data and direct searches suggest that the Higgs boson must be relative light and 
its mass should be roughly in the range of 114.4GeV-208GeV at 95% CL Studying the 
properties of the Higgs potential will reveal details of the mass-generation mechanism in spon- 
taneously broken gauge theories, which can be obtained through measuring the Higgs boson 
self-interactions. Recently, the Higgs boson pair production processes have been widely consid- 
ered, and the cross sections for these processes in the SM have been evaluated at linear colliders 
and hadron colliders. The phenomenology calculation show that it would be extremely difficult 
to measure the Higgs self-coupling Xrhh at the LHC [7], and e^e~ linear colliders, where the 
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study of the e+e~ — ZHH and HHvV can be performed with good accuracy, represent a pos- 
sibly unique opportunity for performing the study of the trihnear Higgs self-couphng P [9l flO]. 
For the center of mass (cm.) energy y/s from 500 GeV up to 1 TeV, the ZHH production with 
intermediate Higgs boson mass is the most promising process among the various Higgs doublet- 
production processes. Since the cross section is relatively large and all the final states can be 
identified without large missing momentum, the process e+e~ ZHH is the best one among 
the various Higgs doublet-production processes to look for the Higgs self-coupling during the 
first stage of the future linear collider. 

We know that the most important Higgs production process at the linear collider is the 
Higgs-strahlung process e^e~ ZH . The correction effects of LH model to this process was 
studied in Ref. [llj. It is found that the correction effects mainly come from the heavy gauge 
boson Bhi in most parameter space, the deviation of the total cross section from its SM value 
is larger than 5%, which may be detected at the future ILC experiment. However, the double 
Higgs-strahlung process e^e~ ZHH includes the trilinear Higgs coupling which is different 
from the process e^e~ ZH . For the process e^e~ —>■ ZHH, the contribution of the LH 
model comes from not only the new heavy gauge bosons Bh, Zh but also the modification of 
the self-couplings of Higgs boson. So the process e^e^ ZHH can also provide some useful 
information about the modification of trilinear Higgs coupling in the LH model to complement 
the study of the process e~^e~ ZH. In this paper, we consider the double Higgs-strahlung 
process e'^e" ZHH and study whether the correction effects of LH model to this process 
can be detected at the future ILC experiment. 

This paper is organized as follows. In section two, we first briefly introduce the LH model, 
and then give the production amplitude of the process. The numerical results and discussions 
are presented in section three. Our conclusions are given in section four. 

2 The process e+e~ ZHH in the LH model 

The LH model is based on the SU (5) / 5*0(5) nonlinear sigma model. At the scale ~ 47r/, 
the global SU{b) symmetry is broken into its subgroup 5*0(5) via a vacuum condensate /, re- 
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suiting in 14 Goldstone bosons. The effective field theory of these Goldstone bosons is parame- 
terized by a non-linear a model with gauged symmetry [SU{2) x f/(l)]^, spontaneously broken 
down to its diagonal subgroup SU{2) x f/(l) which is identified as the SM electroweak gauge 
group. Four of these Goldstone bosons are eaten by the broken gauge generators, leaving 10 
states that transform under the SM gauge group as a doublet H and a triplet $. This breaking 
scenario also gives rise to four massive gauge bosons Bh, Zh and Wh, which might produce 
characteristic signatures in the present and future high energy collider experiments [T^fT^ fTl]. 

After EWSB, the final mass eigenstates are obtained via the mixing between the heavy and 
light gauge bosons. They include the light (SM-like) bosons Zl,Al and observed at exper- 
iments, and new heavy bosons Zh,Bh and Wh that could be observed at future experiments. 
The masses of neutral gauge bosons are given to 0{v'^/p) by 



Ml, = 0, (1) 
Ml = (Ml^Y{l-^^[l + l{c^-sr + l{c'^-s'r-Y]}^ (2) 

f2 2 ( 2 _ 2\2 I 12 2 ^ 2 12 

f2 2 T.(„l2_ I2\2 I 12 2 ^ 2 12 

with X = xh = \gg' ^5g2^^/v2-g'%c2'* ' where v=2AQ GeV is the elecroweak scale, v' is the 
VEV of the scalar SU{2)l triplet and sw{cw) represents the sine(cosine) of the weak mixing 



angle. c{s = yl — cP) is the mixing parameter between SU{2)i and SU{2)2 gauge bosons and 
the mixing parameter c'{s' = — c'^) comes from the mixing between f/(l)i and f/(l)2 gauge 
bosons. Using these mixing parameters, we can represent the SM gauge coupling constants as 
g = gis = g2C and g' = gis' = g2c' . The mass of neutral scalar boson M$o can be given as [12] 



= 4;^. (5) 

The above equation about the mass of $ requires a constraint of 0<x<l, which shows the 
relation between the scale / and the VEV of the Higgs field doublet and triplet (t;, f'). 

Taking account of the gauge invariance of the Yukawa couplings and the U{1) anomaly 
cancellation, one can write the couplings of the neutral gauge bosons Vi{Vi = Z^, Bh, Zh) to 
electrons pair in the form of A^'^*^ = i^^{gy'^^ + g^'^^l^) with [12] 
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The couplings of the gauge bosons to Higgs boson and self-Higgs couphng can be written as 
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In the LH model, the heavy triple Higgs boson exchange can also contribute to the 
process e"'"e^ ZHH . However, compared to the contributions coming from the new gauge 
bosons, the contribution of $° exchanging is suppressed by the order v'^ / which can be seen 
from the couplings between gauge bosons and scalars[T2]. Thus, we can ignore the contribution 
of the scalar triplets to the process e^e~ ZHH . 

The relevant Feynman diagrams for the process e^e~ —>■ ZHH in the LH model are shown 
in Fig.l at the tree- level. 

The invariant production amplitudes of the process can be written as 
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Figure 1: The Feynman diagrams of the process e^e — > ZHH in the LH model. 



Here, G^^(j), M) = ^^^flp is the propagator of the particle. We can see that this process in the 
LH model receives additional contributions from the heavy gauge bosons Zh, Bh- Furthermore, 
the modification of the relations among the SM parameters and the precision electroweak input 
parameters, the correction terms to the SM HHH coupling can also produce corrections to this 
process. In our numerical calculation, we will also take into account these correction effects. 
The main decay modes of Bh and Zh are Vi f f{f represents any quarks and leptons in the 
SM) and Vi — > ZH . The decay widths of these modes have been explicitly given in references 

PI [IE]. 

With above production amplitudes, we can obtain the production cross section directly. 
In the calculation of the cross section, instead of calculating the square of the amplitudes 
analytically, we calculate the amplitudes numerically by using the method of the references [17] 
which can greatly simplify our calculation. 
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3 The numerical results and discussions 



In the numerical calculation, we take the input parameters as M§^ = 91.187 GeV, s^r = 
0.2315 [IH]. For the light Higgs boson H, in this paper, we only take the illustrative value 
Mh = 120 GeV. In this case, the possible decay modes of H are bb, cc, //[l=r, /i or e], gg and 
77. However, the total decay width Th is dominated by the decay channel H bb. In the 
LH model, Th is modified from that in the SM by the order of j p and has been studied in 
Ref. [19]. The cm. energy of the ILC is assumed as y^=500 GeV. 

The absence of custodial SU{2) global symmetry in the LH model yields weak isospin 
violating contributions to the electroweak precision observables. In the early study, global fits 
to the experimental data put rather severe constraints in the / > 4 TeV at 95% C.L.pUl fIT\ . 
However, their analysis is based on a simple assumption that the SM fermions are charged only 
under f/(l)i. If the SM fermions are charged under U{l)i x U{1)2, the constraints become 
relaxed: the substantial parameter space allows / = 1 ~ 2 TeV [22] • If only the U{1)y is 
gauged, the experimental constraints are looser [221 [23] • Therefore , the new contributions 
are suppressed: / = 1 ~ 2 TeV allowed for the mixing parameters c and c in the ranges of 
~ 0.5,0.62 ~ 0.73 [221 [21]- The parameter x < 1 parameterizes the ratio of the triplet and 
doublet VEV's. Taking into account the constraints on /, c, c', x, we take them as the free 
parameters in our numerical calculation. The numerical results are summarized in Figs. (2-4). 

The relative correction 5a /a^^^ is plotted in Fig. 2 as a function of the mixing parameter c for 
f=l TeV, X = 0.5, Mh = 120 GeV, c' = 0.63, 0.67, 0.71, respectively. In Fig.2, Sa = a^ot_^SM 
and a'^" is the tree-level cross section of the ZHH production predicted by the SM. From 
Fig.2, we can see that the absolute value of the relative correction decreases with the mixing 
parameter c increasing. For x = 0.5, the absolute value of da /a^^^ is in the range of 8% — 14% 
in the most parameter space limited by the electroweak precision data. The curves also show 
that with an increase of the value of c', the effect of the LH model is getting stronger. For / < 3 
TeV, the mass of Bh may be lighter than 500 GeV[ll]. In most parameter spaces of the LH 
model, the mass of the heavy gauge boson Zh is larger than 1 TeV. So, there is no s-channel 
resonance effects in our numerical results. 

To see the dependence of the relative correction on the parameter x, in Fig.3, we plot 
5a /a^^^ as a function of the mixing parameter x for c=0.3, d = 0.67 and three values of the 
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Figure 2: The relative correction ba ja^^ as a function of the mixing parameter c for f=l TeV, x=0.5, 
Mil = 120 GeV and different values of the mixing parameter c'. 

scale parameter /. From Fig. 3 we can see that, the absolute value of the relative correction 
decreases as / increasing. The curves also demonstrate that the effect of the LH model is not 
sensitive to x in the range of x < 0.75. This is because the deviations of the cross section 
from the SM are mainly aroused by the contributions of the new gauge bosons when x < 0.75. 




Figure 3: The relative correction da/a^^^ as a function of the mixing parameter x for c=0.3, c'=0.67, 
and / = 1, 1.5, 2 TeV, respectively. 
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However, the figure shows that the absolute values of 5a /a^^^ raised quickly when we take the 
X — > 1 limit and in this case the main contribution to the cross section comes from the Feynman 
diagram involving the trilinear interaction of the SM Higgs boson, which is consistent with the 
conclusions for the contributions of the LH model to Higgs boson pair production at hadron 
colliders [25]. So, if x is large enough the significant correction of the LH model to the trilinear 
Higgs coupling should be observable. The process e^e~ — ZHH can open a unique window 
to probe the Higgs self-coupling which can complement the process e^e" ZH. 
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Figure 4: The relative correction 6a /a^^ as a function of the the scale parameter / for c=0.3, 
c' = 0.67, Mh = 120 GeV, and three values of the mixing parameter x. 

In general, the contributions of the LH model to the observables are dependent on the factor 
1//^. In order to obtain the generic conclusion, we also plot 5a /a^^^ as a function of /(1-3 
TeV) for three values of the parameter x and take c = 0.3, c = 0.68 in Fig 4. One can see 
that the absolute relative correction drops sharply with / increasing, which is consistent with 
the conclusions for the corrections of the LH model to other observables. On the other hand, 
we can see that the absolute relative correction increases as the parameter x increasing. For 
example, the absolute relative correction may reach about 20% when x = 0.8 and / = 1 TeV. 

As has been mentioned above, the total cross section of e~^e~ —>■ ZHH can reach the order 
of 10~^ fb at the ILC. This cross section amounts to about 100 events with the integrated 
luminosity of 1000 fb~^. The la statistical error corresponds to about 10% precision. The 
reference [26] have reviewed the expected experimental precision with which the ZHH cross 
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section can be measured. Even we consider the systemic error of the ILC, the ILC can measure 
the cross section with the precision of 17% assuming a 120 GeV Higgs and the integrated 
luminosity 1000 fb~^ at 500 GeV. The sensitivity can be further improved when a multi- 
variable selection based on a neural network is applied which can reduce the uncertainty from 
17% to 13%. The relative correction of the LH model to the cross section is only comparable to 
the ILC measurement precision and might be detected at the ILC in the favorable parameter 
spaces(for example, small value of /) preferred by the electroweak precision. The statistical 
acuracy to measure the trilinear Higgs coupling is 22% for = 120 GeV with an integrated 
luminosity of 1000/6"^, using the neural network select ion [26]. Only for small / and large x, 
the correction of the LH model to the trilinear Higgs coupling can be detected. 

4 Conclusion 

The little Higgs model, which can solve the hierarchy problem, is a promising alternative 
model of new physics beyond the standard model. Among various little Higgs models, the 
littlest Higgs(LH) model is one of the simplest and phenomenologically viable models. The 
distinguishing feature of this model is the existence of the new scalars, the new gauge bosons, 
and the vector-like top quark. These new particles contribute to the experimental observables 
which could provide some clues of the existence of the LH model. In this paper, we study the 
potential to detect the contribution of the LH model via the process e~^e~ ZHH at the 
future ILC experiments. 

In the parameter spaces(/ = 1 ~ 2 TeV, c = ~ 0.5, c' = 0.62 ~ 0.73) limited by the 
electroweak precision data, we calculate the cross section correction of the LH model to the 
process e^e^ ZHH . We find that the correction is significant even when we consider the 
constraint of electroweak precision data on the parameters. The relative correction varies from 
a few percent to tens of percents. The LH model is a weak interaction theory and it is hard 
to detect its contributions and measure its couplings at the LHC. With the high cm. energy 
and luminosity, the future ILC will open an ideal window to probe into the LH model and 
study its properties. In some favorable case, the relative correction of the LH model to the 
process e^e~ ZHH might be large enough to be measured with high precision at the ILC. 
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Furthermore, the process can also open a unique window into the trihnear Higgs couphng in 
LH modeL 
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